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Abstract 

We introduce a new, approximate sum rule for determining bulk properties of multicom- 
ponent systems, in terms of the pure components properties. We apply this expression 
for the study of lattice parameters, cohesive energies and bulk moduli of binary alloys. 
The correct experimental trends (i.e., departure from average values) are predicted in 
all cases. 
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In this paper we introduce a new expression for the direct calculation of equilibrium 
properties of alloys in terms of the pure components. From this new rule, which can be 
obtained as a particular case from the method for alloys energetics introduced recently by 
Bozzolo, Ferrante and Smith^\ we can simultaneously derive three separate expressions for 
obtaining the lattice parameter, cohesive energy and bulk modulus of alloys. 

Consider an alloy of N different elements A Xa B X]3 ... where x t - denotes the concentra- 
tion of the atomic species i. Let Ec{x), a(x) and B(x) be the cohesive energy, lattice 
parameter and bulk modulus, respectively, of the alloy characterized by the concentrations 
x : (xa, xb, •••)• Then, our approximate expression for the binding energy as a function of 
concentration x and lattice parameter r, J5(x,r), in terms of the corresponding expressions 
for the pure elements is given by 

E(x,r) = Y2 x iEi(r) ( 1 ) 

t 

As mentioned above, this expression, as well as its range of applicability, has a theoretical 
foundation in the recently proposed BFS^ 1 ) method. There is no formal derivation of this 
rule, but, in the framework of BFS, we can find an explanation for its origin and some idea 
about its validity and shortcomings. This will discussed at length elsewhere. Here, our goal 
is to provide three new approximate rules that would enable one to describe equilibrium 
alloy properties in terms of pure materials in a simple way. Although it is approximate, 
recent first-principles calculations closely follow the general behavior predicted by Eq. (1). 
Terakura et al^ 2 ) analyzed the phase stability of several binary alloys systems composed 
of Ni, Pd or Pt as one element and Cu, Ag and Au as the other. From their LDA band 
calculations they found that their predictions for the excess energy A E m (r) for the ordered 
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alloy structures A m B^_ rn {m = 0, ...,4) admit a simple parametrization of the form 


+r. 


( 2 ) 


where the parameters p m ,g m and r m are listed in Ref. 2, with n = 3.5 and m = 0,4 
represent the pure elements.. The excess energy for the disordered alloy can be represented 
by a similar functional form^ 3 ), 


A£ c (r,*)= (y-) 2 ”- (v)" + "' (3) 

with P% n _ E 4 i=oCm(x)p 2n < gn _ c m (x)g^ and R x = Em=o c m(x)r m , where 

c m (x) = x m (l — x) 4-m . If these results behaved according to Eq. (1), then a di- 
rect consequence of this rule, 


ni 

AE D (r,x = — ) = AE m (r) ; m = 0, ...,4 (4) 

would also have to be satisfied. In terms of the coefficients used in Ref. 2. this condition 
can be written as set of simultaneous requirements on these coefficients: 


Pm/ 4 - Pm » Qm/ 4 ~ Qm > ^m/4 “ r m *, m — 1,2,3. (5) 

We show the results of this comparison in Table 1 for the Cu-Ag, Cu-Au and Ag-Au alloys, 
where agreement in the two sides of Eq. (5) is a test of the sum rule. Although the 
comparison is not exact (it’s not supposed to be), the agreement is surprisingly good, 
making the use of Eq. (1) a viable alternative to expressions^ 4 ) for obtaining alloy cohesive 
energies, bulk moduli and lattice parameters, such as Vegard’s law( 5 ) for lattice parameters. 

The application of Eq. (1) is greatly simplified by the fact that the binding energy 
curve of a pure element i, E{(r), can be represented by a Rydberg function^ 6 ) 

E t (r) = -E Cy t(l + a*)e~ a * (6) 
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with a * = q t (r — a,)//,, where a, is the equilibrium lattice parameter of species i and the 
scaling length /,* is related to the cohesive energy of the pure element Ec,i and its bulk 
modulus B{ as if — Ec,i/12nqiBi. qi is the ratio between the equilibrium Wigner-Seitz 
radius and equilibrium lattice parameter of species i. 

Although it is a well-known fact that the binding energy curve (including alloys) is 
of a universal nature^ 6,7 ), the functional form chosen to represent it is just one of many 
choices. If we used this functional form in Eq. (1), a representation of the alloys binding 
energy curve could only be achieved by a numerical fit to a Rydberg function. To avoid this 
complication, and in order to provide simple analytical expressions for approximate values 
of the static properties of the alloys, we expand the Rydberg function to second order on 
both sides of Eq. (1) thus obtaining a simple system of equations for Ec(x), B(x) and a(x). 
The ensuing expressions for these quantities are 


a(x) 


< Qa 2 B > 

< QaB > ' 


B(x) 


1 < QaB > 2 

07< Qa 2 fl >’ 


and 


E(x ) =< E > - < Qa 3 B > + 


< Qa?_B > 

< QaB > 


(7) 

( 8 ) 

(9) 


where 


< 0 >= 5 >, 0 , ( 10 ) 

I 

and Qi = I2n qf. 

If, for example, we focus our attention on binary alloys of elements with the same 
structure as the pure crystals, the Q's cancel in Eqs.(7)-(9). Moreover, it is a simple 
exercise to examine these expressions for their departure from the corresponding average 
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values A 0(x) =< O x > - < 0 > (for 0 = a,B and Ec )• The derivative at the origin of 
A 0 for each one of these quantities are 


/ dAa(x) N 


(a A B A - ubBb ) 


)\ , J 

^dA5(x)\ = _^ (a/i _ afi) 2 

/ x ^=0 < 


and 


V dx 


( dAEc{x) 
\ dx 


x A =0 


= —a A B A (a A — a #) 2 


( 11 ) 

( 12 ) 

(13) 


Similar expressions can be obtained for the limit x A = 1 . Eq.(ll) allows for deviations 
above and below the average values, whereas Eqs. ( 12 ) and (13) predict cohesive energies 
and bulk moduli below the average values for all concentrations. For this simple case, the 
maximum deviation of the lattice parameter from the average value < a > is given by 

67A 


A a r 


(1 + 7)(*U7 + e B) 


(14) 


where 6 — a A — a#, A = a A e A - a#e£, e t = a,i? t and 7 2 = e^/e^. This maximum 
deviation occurs at x A = Finally, it is straightforward to show that Vegard’s law can 
be recovered as a particular case when 7 = 1 . 

The predictions of the trends predicted by Eq. ( 11 ) were successfully tested against 34 
binary alloys of fee (Ag, Al, Au, Cu, 7 -Fe, Ir, Ni, Pd, Pt), bcc (Cr, Cs, a-Fe, K, Li, Mo, 
Na, Nb, Rb, Ta, V, W) elements, as well as Si and Ge in their diamond and fee phases for 
which complete experimental data of the solid solutions is available^ 8 ). Similar results are 
found if other representations of the binding energy curves are used in our proposed rule 
(Eq. ( 1 )). 

The trends in Eq. ( 11 ) depend on a coupling between the lattice parameter and com- 
pressibility as might be expected intuitively. This result is a consequence of the fact that 
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in our expression (Eq. (1)) the binding energy of the pure materials determine the behav- 
ior of the alloys. However, the specific functional form depends on the choice adopted for 
representing the binding energy curve (in this case, a simple quadratic expansion of the 
Rydberg function). The simple example studied in this work predicts the correct trends in 
the slopes (Eq. (11)) of the deviation from the average values as well as the trends in the 
relative deviations for different alloys. 

We have presented a new method for estimating the cohesive energy, lattice parameter 
and bulk modulus for alloys from pure element properties. The method has been successfully 
tested by comparisons with ab-initio calculations and trends in experimental data. 
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System 

m 

^m/4 

Pm 

Qm/4 

Qm 

Rm/4 

1*771 

AgAu 3 

i 

8.46 

8.45 

11.37 

11.35 

1.98 

1.96 

AgAu 

2 

8.33 

8.33 

11.07 

11.07 

1.82 

1.82 

Ag 3 Au 

3 

8.21 

8.21 

10.77 

10.77 

1.66 

1.66 

CuAg 3 

1 

7.84 

7.82 

10.11 

10.08 

1.51 

1.50 

CuAg 

2 

7.60 

7.59 

9.77 

9.79 

1.50 

1.52 

Cu 3 Ag 

3 

7.34 

7.32 

9.42 

9.41 

1.46 

1.47 

CuAu 3 

1 

8.26 

8.22 

11.07 

11.01 

1.95 

1.93 

CuAu 

2 

7.91 

7.88 

10.45 

10.45 

1.77 

1.79 

Cu 3 Au 

3 

7.52 

7.47 

9.78 

9.75 

1.59 

1.59 


Table 1: Comparison of the parameters (P,Q,R) and (p,q,r) as defined by Eqs. (2) and (3) 
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